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Cellular Automata 


• Original idea from John von Neumann and Stanislaw Ulam in the 
I940’s. Originated as an application ofthe mathematical abstraction 
of the growth of crystals to the problem of self-replicating system. 


• Based on cell interaction 


• Conway designed a CA called ‘Game of Life’ in the I970s which is 
the most widely known CA. 
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Cellular Automata principie 


* Cells are positioned in a n-dimensional grid. 

* Each cell is in one of a finite number of States. 


• The cell States are updated each discrete timestep. 

• The cell State update depends on the State of itself and its 
neighboring cells and a transition rule. 
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1D binary Cellular Automata 


Rule: 

• All three cells of the neighborhood white -> white 

• All three cells of the neighborhood black -> white 

• Else -> black 

TO: 

TI: 

T2: 

T3: 

T4: 

T5: 

T6: 
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And thus we get a pattem... 
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A different rule: “Wolfram rule 30 CA” 


00011110(111 to 000 from left to right) 
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Inter esting... 



The “Wolfram rule 30 CA” 
pattern is very similar to the 
patterns of some seashells. 
These seashells exhibit a 
natural variant of CA. 
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Natural 


Each pigment cell secretes pigments 
according to the activity of its 
neighbor cells. The cell band leaves 
the colored pattern on the Shell as it 
grows slowly. 
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Only showing the current State 


lt is also possible to only show the current State of 
the Cellular Automaton: 



And so on... 
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1D Binary CA (1) - Notation 


(j* — (2 } the state °f the CA at time t 

a í G {0,1} - the state of cell i at time t 

2^+1 - the state of neighborhood of cell i at time t 

- the neighborhood radius 


s'i e {°,1} 
r 


0 ! {0,1} 27 +1 —> {0,1} _ the transition rule 




5. 





> 

r 



r 



a* 

a 2 


a 4 

a 5 

a 6 

a 7 

a 8 
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1D Binary CA (2) - Numbers 


• 2" different States for C 



possible States for s. 


r)2r+l 

. 2 2 


possible transition rules 


The result is a huge number of possible rules! 


r 

number of possible rules 

i 

256 

2 

4,3'10 9 

3 

3,4 ■ 10 28 
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1D Binary CA (3) - Boundary cells 


Two possibilities: 


• Toroidal arrangement: a n is linked to a. 



*1 

a 2 

a 3 

a 4 

a 5 

a 6 

a? 

a 8 






• Unconnected boundaries:the border of the grid is enclosed by 
‘zero-cells’ 
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1D Binary CA (4) - Tr ansition rule 


The transition rule 0 can be represented by: 

R = {b 1 ,b 2 ,...,b 2 2r+1 } 


With bj being the output of the transition rule for the mapping of the input 
to 2 2r+l -i. 


Example: R = {0,1,1,1,1,1,1,0} 


Neigborhood State 

111 

110 

101 

100 

011 

010 

001 

000 

Output 

0 

1 

1 

1 

1 

1 

1 

0 


Using the decimal notation ofthe bitstring (Wolfram notation): 
01111110= 126 -> Wolfram rule 126 CA 
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More general Cellular Automata 


The ID binary CA can be extended with: 

• Multiple dimensions 

• Multiple possible States per cell 

• Different neighborhoods shapes 

The effect: 

• More complex behaving CA can be made 

• The number of possible rules grows super-exponentially 
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Multidimensional CA transitions 


The transition rules are defined the similar way: 

• The cells of the neighborhood are numbered 

• The transition still rule maps a neighborhood State to a new State: 

0: {0,1 } n -> {0,1}, n = S(d,r) 

• S(d, r) depends on: 

- The number of dimensions 

- The radius 

- The neighborhood shape 

ryS(d,r) 

• The number of possible rules becomes 2 
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Multidimensional CA neighborhoods 


Two main neighborhood shapes: 


• Von Neumann neighborhood 


• Moore neighborhood 


Different neighborhood shapes could be used as well. 
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Von Neumann and Moore radius 


R=1 R=2 R=3 


• Von Neumann 
neighborhood 


Moore 

neighborhood 
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Conway’s Game of Life (1) 


• A 2 dimensional binary CA that uses a Moore neighborhood and 4 
general rules 

• Devised by the British mathematician John Horton Conway in the 
I970s 

Rules: 

Any live cell with fewer than two live neighbours dies of loneliness 
Any live cell with more than three live neighbours dies of overcrowding 
Any live cell with two or three live neighbours lives on 
Any dead cell with exactly three live neighbours comes to life 
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Conway’s Game of Life (2) 
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Conway’s Game of Life (3) 



Strudel 


Pulsator 



Extraordinary cell interaction 


Jellyfish 


All around: There exists quite a few 
patterns in Conway’s Game of Life that 
show remarkably lifelike behavior. 


Octagon 

ÍS 


Swimmer 
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Four classes of behavior 


1 . the CA settles into a stable State. 

2. the CA exhibits simple periodic behavior (oscillation). 

3. the CA acts chaotically (like a turbulent fluid). 

4. the CA exhibits behavior between the classes 2 and 3. It is 
behavior is more complex than simple oscillation but more 
orderly than chaos.Therefore this class of behavior is often called 
the “edge of chãos”. 

Class 4 CA are the most interesting CA and are considered suitable 
for computation because they exhibit long-range interactions (which 
can be used for communication), numerous semi-stable States 
(which can be used for memory), and other desirable characteristics. 
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Wolfram Class 1 


• Develops into a stable State 

• Complete loss of information 
about initial State 


• Corresponds to a dynamical 
system with a single attractor 
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Class 1 State transition 
function 
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Wolfram Class 2 


• Periodical structures 

• Finite-size automata are 
automatically in class I or 2 


• Cyclic, dynamical system 


t 
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Class 2 State transition 
function 
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Wolfram Class 3 


• No cycle of finite lenght 

• Very sensitive to changes of 
initial configuration 


• Deterministic chãos 
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Class 3 State transition 
function 
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Wolfram Class 4 


• Local structures travei and 
interact with each other 

• Information transmissions 
allows for complex 
computations 


• Between cyclic systems and 
deterministic chãos 
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Class 4 State transition 
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Classification Methods? 


• Membership of a particular class: Undecideable 


• Equivalent to halting problem 


• Langdon-Factor gives an indication of class membership, though. 
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Artificial Life in CAs 


• Universe 

• Physical laws 


• Temperature 


• Biomolecules 
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Natural 


* -CA 

* • Transition function 

* • Langdon-Factor 

> • Periodic patterns 
(virtual automata) 
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Definition: Langdon-Factor 


• Ground State' * ^ neighbour cells of a cell c are in this State, then 

c itself attains this State at the next time step. 

• Number of neighbouring cells: 

P 

• Possible number of cell States: 

k 

• Number of neighborhood States: 

k p 

• Number of transitions into ground State: 


k p -n 

Definition: À = q 

k p 
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Langdon-Factor and Wolfram Classes 


• Simulation: 
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Reversibility and unreachable States 


• A Cellular Automata is reversible if for every current State of the 
CA there is exactly one past State (preimage) 

• Cellular Automata that are not reversible have patterns for which 
there are no previous States.These patterns are called ‘Garden of 
Eden patterns’ 

• For ID CA there exist algorithms that can find preimages, and 
decide whether a rule is reversible or irreversible 

• For CA of two or more dimensions it has been proved that the 
reversibility is undecidable for arbitrary rules 
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The inverse problem 


• Given a number of pre-selected global properties, find a cellular 
automata rule that will have these properties 


• The inverse problem of deducing the local rules from a given global 
behavior is extremely difficult because of the huge search space (all 
the possible rules) 
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So what? 


• Certain rules tend to exhibit organizational behavior independently 
of the initial State of the CA. 

• The organizational behavior shows a remarkable resemblance to 
many real-life processes and can be used to simulate these real-life 
processes. 

• The organizational behavior suggests that there is some form of 
communication going on between the cells of the CA that goes 
beyond the scope of the cell neighborhoods. 


Simple cells can communicate and work together to form a 
complex system capable of powerful behavior! 
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Simulation with Cellular Automata 


Cellular Automata can be used to simulate: 

• Biological processes 

• Câncer cells growth 

• Forest fires 

• Social movement 

• ... more ... 

Left: the forest fire example shows that CA can be 
used to obtain very realistic simulations. 

Source: 

http://www.eddaardvark.co.uk/svg/forest/forest.html 
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Finding rules / solving problems 


Cellular Automata could be used for solving problems in a parallel 
computational manner. 

Consider the following problems: 

• Majority problem 

• AND problem 

• OR problem 


Goalrfind the rules that give the correct answer in most of the cases 
(preferably 100% correct). 

Because of the huge number of possible rules, this search is far from 
trivial! 
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Solving problems (1) - Majority 


Let: 

X = the relative number of ones in C 0 
Task: 

A. > 0.5 iterate to ‘all ones’ State 

À < 0.5 iterate to ‘all zeros’ State 



O Position 148 O Position 148 

Also possible in multiple dimensions: 2D Majority problem 
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Solving problems (2) - AND Problem 


Given the top-left and bottom-right cell v, and v 2 in a 2D binary CA 
with unconnected boundaries: 


if (V| AND v 2 ) -> iterate to ‘all ones’ 
else -> iterate to ‘all zeros’ 

( 0 , 0 ) 

( 0 , 1 ) 

( 1 , 0 ) 

( 1 , 1 ) 
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Solving problems (3) - XOR Problem 


Given the top-left and bottom-right cell v, and v 2 in a 2D binary CA 
with unconnected boundaries: 


( 0 . 0 ) 

(0,1) 

( 1 , 0 ) 

( 1 , 1 ) 


if (v, XOR v 2 ) 
else 


-7 iterate to ‘all ones’ 
-> iterate to ‘all zeros’ 
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Cryptography 


• The inverse problem of deducing the local rules from a given global 
behavior is extremely difficult because of the huge number of 
possible rules. 


• This property gives rise to the question whether CA can be applied 
for cryptographical use in public key cryptography. 

• Wolfram suggested Rule 30 as a possible stream cipher for use in 
cryptography. 
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Evacuation Simulation 


• Goal 

- Optimal sizing and placement of doors 

• A pproach 

- CA simulation of evacuation behaviour of humans 

• Background 

- CAs have been used for evacuation simulation for a long time 

- Consideration of obstacles 

- Irrational behaviour of humans? 

“Cellular Automaton Model for Evacuation Process with Obstacles”, 

A. Varas et al., 2007 
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CA Description 


• 2-D, rectangle = room to be evacuated 

• 3 possible States: 

- Empty 

- Occupied by person 

- Occupied by obstacle 

• Transition function: 

- Everybody wants to go to the exit - shortest path 

- If two persons want to enter the same cell, a random 
decision is made 

- Persons can fali down / not move anymore 


Prof.Thomas Bàck 


Natural Computing Group 


Cellular Automata 39 



No Obstacle 































































o 





















o 

o 


















o 

o 

o 

o 

















o 

o 

o 

o 
















O 


o 

o 

o 

o 
















o 


o 


















o 

o 

o 

o 

















o 

o 



















o 




































































































Prof.Thomas Bàck 


Natural Computing Group 


Cellular Automata 40 















































Obstacle 
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CA for Urban Growth Modeling 


• Goal 

- Find transition rules for simulating city growth behaviour 

• A pproach 

- Transition rules are parameterized templates 

- Best parameter values are found using genetic algorithms 

- This is called calibration (to real-world data) 

• Background 

- CA is used for modeling growth of cities 

- How to find optimal transition rules? 

- Use parameterised rules and adjust parameters by using 
evolutionary computation / optimization 

“Genetic Algorithms for the Calibration of Cellular Automata Urban 
Growth Modeling”, J. Shan et al., 2008 
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Approach 



CA state 1987 


Parameter values for 
Transition rules 


O 



Initial state 


o 



Fitness of parameter 
values 
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1 step of simulation = 1 year 
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Compare 
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Further reading 


General Cellular Automata: 

• http://en.wikipedia.org/wiki/Cellular_automata 

• http://cell-auto.com/ 

Conway’s Game of Life: 

• http://www.math.com/students/wonders/life/life.html 

• http://en.wikipedia.org/wiki/Conway's Game of Life 

A 3D version of Conway’s Game of Life: 

• http://fbim.fh-regensburg.de/~saj39 122/doefe/ 
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Questions? 



Universiteit 

Leiden 

The Netherlands 
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